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Q ■ Abstract 
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^ . This paper surveys and unifies a number of recent contributions that have collectively developed a metric for 

decentralized wireless network analysis known as transmission capacity. Although it is notoriously difficult to 
' derive general end-to-end capacity results for multi-terminal or ad hoc networks, the transmission capacity (TC) 

framework allows for quantification of achievable single-hop rates by focusing on a simplified physic al/MAC- 



layer model. By using stochastic geometry to quantify the multi-user interference in the network, the relationship 
between the optimal spatial density and success probability of transmissions in the network can be determined, 
and expressed - often fairly simply - in terms of the key network parameters. The basic model and analytical 
tools are first discussed and applied to a simple network with path loss only and we present tight upper and lower 
bounds on transmission capacity (via lower and upper bounds on outage probabihty). We then introduce random 
channels (fading/shadowing) and give TC and outage approximations for an arbitrary channel distribution, as 

o 

' well as exact results for the special cases of Rayleigh and Nakagami fading. We then apply these results to show 

how TC can be used to better understand scheduling, power control, and the deployment of multiple antennas 
OO . in a decentralized network. The paper closes by discussing shortcomings in the model as well as future research 
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directions. 



I. INTRODUCTION 

This paper presents the recently developed framework for the outage probability and transmission 
capacity [[T| in a one hop wireless ad hoc network. The transmission capacity is defined as the number 
of successful transmissions taking place in the network per unit area, subject to a constraint on outage 
probability. In addition to being of general interest, the advantange of transmission capacity - relative 
to, say, the transport capacity or average sum throughput - lies largely in that it can be exactly derived in 
some important cases, and tightly bounded in many others, as we shall show. From the expressions and 
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approach given in this paper the exact dependence between system performance (transmission capacity, 
outage probability) and the possible design choices and network parameters are laid bare. In contrast 
to the proposed framework, nearly all other work on ad hoc networks must resort to scaling laws or 
numerical simulations, in which case intuition and/or precision is usually lost. 

The first goal of this paper is to concisely summarize the new analytical tools (largely drawn from 
the field of stochastic geometry (21, lEl) that have been developed over numerous papers by the authors 
and others. Because these techniques have been developed somewhat independently depending on the 
problem of interest, the system model in ^|ll]applied to the baseline model of pathloss attenuation without 
fading in will help newcomers to the area understand the various approaches in context. 

The second goal is to show how this framework can be used to give crisp insights into wireless 
network design problems. In the past few years, the transmission capacity approach has been applied 
to various design problems by a growing group of researchers (see [[II, [|4l-[[8l). Although transmission 
capacity was originally developed to analyze spread spectrum in ad hoc networks, it has proven to be a 
metric with considerable breadth of application. Since decentralized wireless networks are generally very 
difficult to characterize, the intuitive and simple-to-compute qualities of transmission capacity have made 
it a popular choice for a large number of possible systems, including: i) direct-sequence and frequency- 
hopping spread spectrum [[fl, [SI, [[U, ii) interference cancellation [[H, [[TOl . in) spectrum sharing in 
unlicensed, overlaid, and cognitive radio networks [[6l, [[3, [[TT|. [[T2|. iv) scheduling [[TOll and power 
control [[T3]| . [[T4l . v) and the use of multiple antennas (which had resisted characterization by other 
methods) [[H, [[T5l - [|2TI . Other researchers have also further studied the basic tradeoffs between outage 
probability, data rate, and transmission capacity for general networks [[22|. We selectively discuss some 
of these applications. JV| addresses networks with fading channels, with a focus on Rayleigh ( §IV-BI) 
and Nakagami ( §IV-CI) fading, scheduling ( §IV-DI) . and power control ( §IV-EI) . fVl addresses the use 
of multiple antennas, with discussions of diversity ( §V-AI) . spatial interference cancellation ( §V-BI) . and 
spatial multiplexing ( §V-CI) . 

The third goal of the paper is to stimulate new efforts to further the tools presented here, both in 
making them more general and in applying them to new problems. We readily concede that the presented 
model has some nontrivial shortcomings at present, and we identify those as well as possible avenues 
forward in fVT[ 

II. System model 

We introduce the system model in §II-A[ discuss relevant mathematical background in §II-B[ and 
elaborate on the connection with transport capacity in §II-C[ 
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A. Mathematical model and assumptions 

We consider an ad hoc wireless network consisting of a large (infinite) number of nodes spread 
over a large (infinite) area. The network is uncoordinated, meaning transmitters do not coordinate with 
each other in making transmission decisions. That is, nodes employ Aloha |[23ll {i.e., in each slot, each 
node independently decides whether to transmit or to listen) as the medium access control (MAC) 
protocol. We view the network at a snapshot in time, where the locations of the transmitting nodes at 
that snapshot are assumed to form a stationary Poisson point process (PPP) on the plane of intensity 
A, denoted n(A) = {Xi], where each Xi G is the location of interfering transmitter i. The PPP 
assumption for node locations is valid when the uncoordinated transmitting nodes are independently and 
uniformly distributed over the network arena, which is often reasonable for networks with indiscriminate 
node placement or substantial mobility. If intelligent transmission scheduling is performed, the resulting 
transmitter locations will most certainly not form a PPP, so this paper's analytical framework is primarily 
applicable to uncoordinated transmitters. Although suboptimal, such a model may be reasonable in cases 
where the overhead associated with scheduling is prohibitively high, for example due to highly mobile 
nodes, bursty traffic, or rigid delay constraints. We also note that this framework has been extended to 
CSMA, and the gains are not that large over Aloha [|24l . [|25l . Viewing the network at a single snapshot 
in time restricts our focus to characterizing the performance of one-hop transmissions with specified 
destinations. That is, our attention is on (uncoordinated) MAC layer performance, but our model neither 
addresses nor precludes any multi-hop routing scheme. These model limitations are further discussed 
in gVll 

Each transmitter is assumed to have an assigned receiver at a fixed distance r (meters) away. 
This assumption may be easily relaxed {e.g., see [[T3l and lfT4l ) but at the cost of complicating the 
derived expressions without providing additional insight. The set of receivers is disjoint with the set of 
transmitters. Because the network is infinitely large and spatially homogeneous, the statistics of n(A) 
are unaffected by the addition of a placed transmitter and receiver pair, and, more importantly, this pair 
is "typical" in that the performance experienced at the reference pair characterizes the node-average 
performance in the network (Slivnyak's Theorem f2\). Without loss in generality we place the reference 
receiver at the origin (o), and the reference transmitter is located r meters away. See Fig. [TJ Note that 
the locations of the other receivers are not important because the reference receiver's performance only 
depends upon the positions of the transmitters. 

Each transmitter is usually assumed to employ unit transmission power (except when we discuss 
power control in §IV-EI) . The channel strength is assumed to be solely determined by pathloss and 
fading, i.e., the received power at distance d is Hd~°', where a > 2 is the pathloss exponent and H is 
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the fading coefficient. All fading coefficients are assumed to be independent and identically distributed 
(iid). This simplified model has been shown to capture the key distance dependency in ad hoc networks, 
and minor alterations to it such as adding an attenuation constant or forcing the received power to be 
less than one increase the analytical complexity with little apparent benefit [26] . We study networks 
without fading {H = 1) in ^ fllll then with fading in ^ ffVl 

We treat interference as noise, assume that the ambient/thermal noise is negligible, and assert trans- 
mission success to be determined by the signal to interference plus noise ratio (SINR) lying above a 
specified threshold (3. The assumption of negligible thermal noise may be easily relaxed (e.g., see lfT3l 
and [[Ml) but at the cost of complicating the derived expressions without providing additional insight. 
The outage probability (OP), denoted by q, is the probability that the signal to interference ratio (SIR) 
at the reference receiver is below a specified threshold (3 required for successful reception: 

,(A) . P(siR < « ^ P { j:J;^^x.i~^ > • « 

where Y = Z]ien(A) -^il^d ° defined as the aggregate interference power seen at the reference 
receiver at the origin, normalized by the signal power Sr^". The last expression in ([T]) highlights the 
fact that, conditioned on S, the OP is the tail probability of the aggregate interference level expressed 
as a shot noise process. 

The randomness is in the interferer locations, {Xi}, and the fading coefficients, S and {h}. The OP 
is a function of a,(3,X,r and the fading statistics. Note that q is continuous monotone increasing in 
A and is onto [0, 1]. Our primary performance metric is the transmission capacity (TC) which takes a 
target OP e as a parameter: 

c{e) = q-\e){l-e), eG(0,l). (2) 

It is the spatial intensity of attempted transmissions q^^{e) associated with OP e, thinned by the 
probability of success, 1 — e. The quantity e is a network-wide quality of service measure, ensuring a 
typical attempted transmission will succeed with probability 1 — e. The transmission capacity has units of 
number of transmission attempts per unit area, i.e., it is a measure of spatial intensity of transmissions. 
Note that the OP q{\) is defined for an arbitrary transmission intensity A, and c(e) is simply that value 
of A such that g(c(e)/(l — e)) = e. The definition of TC is motivated by several factors: i) fixing the 
OP at g = e is a useful and simple, albeit coarse, characterization of network performance, ii) the TC 
is tractable and can be computed, or at least bounded, for many useful network design questions. A 
summary of the mathematical notation employed in this paper is given in Table [H 
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B. Mathematical background 

The key underlying mathematical concept is the shot-noise process first developed in 1918 [|27ll . 



where {tj} is a stationary Poisson point process (PPP) on M and h{t) is a (linear, time-invariant) impulse 
response function ll28l . Il29l . Here Y{t) is the superposition of all signals, appropriately attenuated to 
time t. If we instead interpret {tj} as locations on the plane, t as the location of a reference receiver, h{t) 
as a channel attenuation function, and t — tj as the distance from tj to t, then Y{t) may be interpreted 
as the cumulative interference power seen at t. A power-law impulse response, h{t) = Kt~" ll30l makes 
the process {Y{t)} Levy stable ll3T1l . 

The use of spatial models in wireless communications dates back to the late 1970's [|32l . [|33l . There 
was in fact quite extensive work on the model in which nodes are located according to a 2-D PPP, Aloha 
is used, a routing protocol determines the node for which each transmitted packet is intended for, and 
the received SINR and specifics of the communication protocol determine conditions for transmission 
success; see ll34l for an overview of early results. The aggregate interference process in an ad hoc 
network was first recognized as Levy stable in Il35l - ll37ll . and its characteristic function was studied in 
ll38l . A series of papers by Baccelli et al. demonstrated the power of stochastic geometry for modeling 
a wide range of problems within wireless communications, as summarized in [|39l , ll40l . 

We note that there have been several very helpful tutorials on applying stochastic geometry to wireless 
networks developed in the last year, including the comprehensive two-volume monograph by Baccelli 
||39l , a monograph by Ganti and Haenggi that has many of the available results on non-homogeneous 
Poisson node distributions [(41 J. a summary tutorial article for a JSAC special issue on the topic [|42||. 
and a tutorial by Win et al. on characterizing interference in Poisson fields [|43ll . We refer readers to 
those references (and flU, [[3l) for background. 

C. Relationship to transport capacity 

The general subject of the paper is the analysis of capacity and outage probability of wireless ad hoc 
networks. Ideally, one could determine the capacity region of an ad hoc network, which would be the 
set of maximum rates that could be achieved simultaneously between all possible pairs in the network, 
and hence is n{n — 1) dimensional for n (full-duplex) users. Even if this was obtainable - which it 
has not been despite considerable efforts ll44l - it would still likely not capture some key aspects of an 
ad hoc network, which call for information to be moved over space. Gupta and Kumar pioneered an 
important line of work on transport capacity in [|45l . which measures the end-to-end sum throughput 



oo 




(3) 



6 



of the network multiplied by the end-to-end distance. Representative publications include [|46ll - [l50ll . 
A key feature of all these works is that it is not possible to compute the exact transport capacity in 
terms of the system parameters, and although bounds and closed-form expressions are available in some 
cases, the best-known results are stated in the form of scaling laws that quantify how the volume of 
the capacity region grows with the number of nodes in the network. The most accepted conclusion is 
that the capacity grows sublinearly as Q{^/n), which can be achieved with multi-hop transmission and 
treating multi-user interference as noise, as proven in several different ways P31 . ||49ll , [ISTIl including 
recently using Maxwell's equations [|52ll . Generous assumptions on mobility ll53l . bandwidth [|54ll . or 
cooperation [|55l result in more optimistic scaling laws. 

The transport capacity, Cxin), is defined as the maximum distance- weighted sum rate of commu- 
nication over all pairs of n nodes HSI . In an extensive network, where the density of nodes per unit 
area is constant, the transport capacity has been shown to grow as CT{n) = Q(n) as n oo, with 
units of bit- meters per second [|46ll . Roughly speaking, there can be Q(n) simultaneous nearest-neighbor 
transmissions in the network, and the distance and the rate of communication between nearest-neighbors 
are both 0(1), yielding Crin) = 0(n). 

Comparison of transport capacity and TC is facilitated by normalizing the transport capacity by the 
network area, A(n) = Q{n), giving Cr(n)/y4(n) = 6(1) in units of bit-meters per second per unit 
area. Within the TC framework, assuming communication at the Shannon rate of log2(l + /?), this 
metric also is 6(1) and is precisely c(e)log2(l + (3)r. Thus, transmission and transport capacity are 
consistent in the scaling sense. Furthermore, by abstracting out the end-to-end and multihop aspect of 
the network, the transmission capacity framework allows for a detailed study of the critical constant 
term; this is generally very difficult to do if using transport capacity. Transport capacity and TC are 
complementary metrics: transport capacity gives order optimal throughput, optimized over all MAC and 
routing techniques, while TC gives detailed performance and design insights for the lower layers of the 
network. 

III. Baseline model: Path loss only 

In this section, a baseline model is presented where the only randomness is in the position of the 
nodes, i.e., there is no fading (S = I and h = I for each i in ([T])). Upper and lower bounds are given on 
outage probability and transmission capacity, emphasizing the impact that dominant (strong) interferers 
have on the sum of the interference. The impact of fading is addressed in 
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A. Exact results 

The points of the 2-D PPP of intensity A, i.e., n(A) = {XJ C M?, may be mapped to a 1-D PPP of 
unit intensity using Corollary 2 in [56J. In particular, 7rA|Xjp ~ Tj, where is the squared distance 
from the origin of the i^'^ nearest transmitter, and Tj is the distance from the origin of the i*^ nearest 
point in a unit intensity 1-D PPP. Applying this to the normalized interference power F in ([T]) gives: 

iGn(A) ien{A) jeni(i) 

where the notation ni(l) indicates a 1-D PPP of intensity 1. The corresponding OP in © becomes 



\ i6ni(i) J 



g(A) = P I (vrr^A)t E > ^ I = ^ > (^^^^f^ ) = ^z. (((-^A)t/3)-^j , (5) 

where = X]ieni(i)-^i ^ is a random variable whose distribution depends only on a and Fz^{-) is 
the complementary cumulative distribution function (CCDF) of Z^. Using F^^(-) to denote the inverse, 
and solving Fz^ ^((vrr^A)^/?) ^ j = ^ for ^ allows the TC to be written as: 

= ^ ,J -■ (6) 

These transformations highlight that the essential difficulty in computing the OP and the TC lies in 
computing the distribution of the stable rv Za. 

In fact the only a > 2 for which Za has a distribution expressible in closed-form is for a = 4, which 
is the inverse Gaussian distribution. Important early results for this special case are due to Sousa and 
Silvester [|36l (Eqn. (21)). In particular, they give an exact expression for the OP in terms of the CDF 
of the standard normal rv, Q{z) = F{Z < z), for Z ~ N{0, 1): 

g(A) = 2Q (v^ATrr^y;^) - 1. (7) 

The corresponding exact expression for the TC is: 



_ v^(i^e)g--((i + .)/2)^ 

An additional exact result is given for the case of Rayleigh fading in ^IV-A[ The general unavailability 
of closed form expressions for the distribution of Z^ motivates the search for lower and upper bounds, 
which we discuss next. 
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B. Lower outage bound: dominant nodes 

A lower bound on the probability of outage is obtained by partitioning the set of interferers 11 into 
dominating and non-dominating nodes. A node i is dominating if its interference contribution alone is 
sufficient to cause outage at the receiver. We call dominating nodes near (n) nodes and non-dominating 
nodes far (f) because because dominating nodes must be within some distance of the origin, and non- 
dominating nodes must be far from the origin. The dominating nodes may be defined geometrically as 
the interferers located inside a disk centered at the origin of radius /3^r: 

n"(A) = |x, : < /?! = |X, :|X,|</3^r}=n(A)n6(o,/3^r). (9) 

Here b{o, d) = {x eM? : \\x\\ < d} denotes the ball centered at the origin o of radius d. The aggregate 
interference, normalized by the received signal power r"", may be split into aggregate dominant and 
aggregate non-dominant interference: 



jen(A) jenn(A) j^n"(A) 

where Y = + Y^. The lower bound is obtained by ignoring the non-dominant interference: 

q{X) = F(Y^ + Y'>^^ >f(y->^^ =g'(A). (11) 

Note that, by construction, the event {Y^ > ^} is the same as the event {n"(A) ^ 0}, which is simply 
the complement of a void probability for a Poisson process: 

g'(A) = 1 - P(n'^(A) = 0) = 1 - e^^K-^'^-OI = 1 - ^-^-r^^i _ ^^^2) 
By solving g'(A) = e for A we get an upper bound on q^^{e), which yields a TC upper bound: 

c^ie) = i^-^)^o^i\-^y' = \ ^ + o(e^) as e 0. (13) 

The right hand side is obtained by observing that the first order Taylor series expansion of (1 — e) log(l — 
e)^^ around e = equals e + O(e^), where O(-) is the standard "big-oh" notation ||57l . Neglecting the 
O(e^) term gives an error e — (1 — e) log(l — e)^^ ~ 0.005 for e = 0.1. The right hand side may be 
interpreted as a disk packing statement. In particular, the maximum number of transmissions per square 
meter for fixed a, /3, e, r is found by packing disks of radius R{a, (3, e, r) = each disk with a single 
transmitter at the center. This radius clarifies the dependence of the supportable density of transmissions 
on these four key model parameters. 
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C. Upper outage bounds: Markov, Chebychev, and Chernoff bounds 
We may decompose the outage event in (fTTI) as: 

,(A) = p ({y- > i} u {y > i} u {y- < i.r' < i,y + y' > i}) . (14) 

In words: the event {Y^ + Y^ > 1/(3} means either or Y^ individually exceed 1/(3, or they are both 
below 1/(3 but their sum exceeds 1/(3. By construction, however, the event [Y^ < 1/(3} is the same as 
the event {F" = 0}, which means the third event in (O is null. The probability of the remaining first 
two events may be written as: 



q(X) = P > 



(3 



p ( > 



(3 



-P 



Y->^y (y' > = q\X) + (1 - q\X))F (y' > , 



(15) 

where we have exploited the independence of Y^,Y^ and applied the definition of g'(A) in (fTTI) . 
Substituting (fT2l) for g'(A) into (fTSi) . we obtain an upper bound on g(A) by an upper bound on 
P (Y^ > 1/(3). We presently give three such bounds, using the Markov and Chebychev inequalities 
and the Chernoff bound. Although the details of the analysis below differ for each of the three bounds, 
the general techniques is the same: upper bound P (Y^ > 1/(3) using the inequality, substitute into (fTSi) . 
then seek a simple expression that upper bounds the resulting expression. 

The Markov inequality [58 1 gives F(Y^ > 1/(3) < (3E[Y^]. Campbell's Theorem O states that if 
{Xi} are points drawn from a PPP of possibly varying intensity X{x) then 



E 



ien 



f{x)X{dx). 



(16) 



Applying this to find E[F^] is straightforward after a change of variable to polar coordinates: 

1 



E 



t-"A27rMt = A = uA, 

a — 2 



(17) 



jennb(o,s) 

where s = (3^r. Multiplying (fTTI) by /3 and combining with (fTSi) . an upper bound on outage is 



q{X) < q 



,u, Markov 



(A) 



-Xnr'^f} a 



-A. 



a 



(18) 



Using the bounds 1 — e ^ < A and e ^ < 1 for A > and simplifying gives a "relaxed Markov" upper 
bound: 



a 



-Txr^fi^X. 



a-2 a-2 

Setting (fT9l) equal to e and solving for A gives a relaxed Markov lower bound on the TC: 



^i,Markov/ \ 



a-2 



a 



Tir'^(3- 



+ O(e^) as e ^ 0, 



(19) 



(20) 
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which is clearly smaller than the TC upper bound of (fT3l) by a factor (a — 2) /a. The right hand side 
is obtained by observing that the first order Taylor series expansion of e(l — e) around e = equals 
e + O(e^). Neglecting the O(e^) term gives an error e — e(l — e) = = 0.01 for e = 0.1. 
Campbell's Theorem also gives the variance of the far-field aggregate interference: 



Var(rO = E 



r 

ieunb{o,s) 



Ar^" r r'^27ftdt = '''''^"' a = a'X (21) 

a — 1 



We use (flTl) and (|2T)) and Chebychev's inequality [58] on the far-field aggregate interference (assuming 

E[Y'] < i), as: 



f(y^ <¥(^\Y^ - E[Y^] I > ^ - Efyf]^ < — 



^ (22) 

^ -/iA 

Substituting (|22)) into (fT5l) and using the bounds 1 — e'"^ < A and < 1 for A > and simplifying 
gives a "relaxed Chebychev" upper bound: 



^«,Chebychev^^^ < TTT^ X + ^^^^ ^. (23) 



a-2 

This expression is quadratic in A; setting equal to e and solving for A gives the relaxed Chebychev 
lower bound on the TC. 

The Chernoff bound [l58ll may be used to obtain an upper bound on the OP: 

P (y^ > ^) - j^J^ = ^"^P (^^ ~ /i (e^""^"" - l) a;da;^ | . (24) 

This expression may be obtained by computing the moment generating function of Y^ restricted to 
b{o,v) and then letting f — )■ cxo, as in IfTOl . [|36ll . The final upper bound on OP is then: 

g"'C"(A) ^ 1 - (^1 - exp |- sup Q - 27rA (e^''"^"" - l) xdx^ }^ ^-x^r^^^ _ (25) 

Although the Chernoff OP upper bound is in some cases tighter than its Markov or Chebychev coun- 
terparts, it depends upon A in a complicated way which precludes a closed-form expression for the 
corresponding lower bound on the TC. In this case, numerical inversion techniques must be applied. 
Sample lower and upper bounds and exact expressions for both OP and TC are shown in Fig. |2l 



D. Tightness of the lower bound: sub-exponential distributions 

Comparing the lower outage bound (fT2l) with the upper outage bound (|25l) . and glancing at Fig. [2l it 
is apparent that the (simple) lower outage bound is much tighter than the (complicated) upper bound. 
One explanation for this comes from the fact that the random interference contribution of each node 
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obeys a subexponential distribution \59l . Consider n points distributed independently and uniformly 
over a disk of radius d centered at the origin, denoted {Xi, . . . It is straightforward to establish 

the CCDF of the individual interference rvs, V = to be Fv{v) = (^vid^ for v > d""^. 

A sufficient condition for a distribution to be subexponential is that limsnp^^^vhyiv) < oo where 
hv{v) = (— logFv(f)) is the hazard rate function. In our case, we find vhv{v) = ^, ensuring Fy 
is subexponential. A defining characteristic of subexponential distributions is the fact that sums of iid 
rvs {Vi, . . . , Vn} typically achieve large values v by having one or more large summands (as opposed 
to a large number of moderate sized summands) ll59l : 



lim — . . =1, n>2. (26) 

maxjKi, ...,Vn}>v) 



V—^CX3 



Because the interference contributions from each node are subexponential, it follows that the probability 
of an outage event {Vi + ■ ■ ■ + Vn > v} (for large v) approximately equals the probability of there being 
one or more dominant nodes with Vi > v. Replacing J2ieu{\) l^d " O with gives 
y = r~"^. Thus V is large if either (5 is small (receiver can decode small SIR) or r is small (Tx and 
Rx are close together) . For small v (meaning both (3 and r are large), outage occurs more easily, and 
in particular, outage may occur due to the aggregate interference being large, even though there may 
not be any dominant nodes. This argument holds for fixed d and n, but gives intuition as to why the 
dominant interference lower bound is tight. 

E. Optimization of SINR Threshold and Outage Constraint 

The SINR threshold (3 and the outage constraint e, which are treated as constants in the TC framework, 
are generally under the control of the system designer and should be chosen reasonably. A meaningful 
objective is maximization of the area spectral efficiency c(e) log2(l + (3), i.e., the product of successful 
density and spectral efficiency. Using Q, the joint maximization over e) can be written as: 

(Fz\e))~^ (1 -e) 

max c(e) log^ll + /3) = max ^ ^" ' log2(l + /?). (27) 

This clearly allows for separate maximizations of (3 and e: 

r = argmax i^^^^^^, = argmax (F^J(e))-^ (1 - e), (28) 

where the optimizers (3* and e* depend only on the path-loss coefficient a. In |[TT1 Section IV], where 
a related but slightly different problem is studied, a closed-form solution for (3* was found: 

/3* = et+M-f^-^)-l (29) 
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where yV{z) is the principle branch of the Lambert W function. Fz^{-) is not known in closed form, 
and thus e* must be determined numerically. In Fig. [3l (3* and e* are plotted versus a, and both are 
seen to be increasing in a. 13* is consistent with normal operating spectral efficiencies, while e* shows 
that the optimal e that maximizes the TC may be unacceptably large. Although such a large outage 
provides a large area spectral efficiency, it also translates directly to long transmission delays and energy 
inefficiency. This analysis highlights a key drawback in unrestricted (spatial) throughput maximization: 
the max-throughput operating point may have an unacceptably high associated OP. The TC framework 
captures this tradeoff by definition: it gives the maximum spatial throughput subject to a specified OP 
constraint. 



We now evolve the discussion to consider channels that also have a random fluctuation about the 
path loss, commonly known as fading or shadowing. The SIR in ([U) models the scenarios discussed in 
this section where random variable S represents the desired signal fade and Jj the fading coefficient 
from the i-th interferer. We assume S is drawn according to some distribution Fs and each Jj according 
to Fi with S, Ii, I2, . . . independent. Independent fading is assumed for tractability; computing the OP 
and TC in correlated fading will be more difficult. 

We first develop a framework for analyzing OP and TC with an arbitrary random channel, and then 
show exact results on OP and TC for Rayleigh and Nakagami fading. It is initially surprising that exact 
results on OP and TC can be computed with certain types fading, but not without fading; recall in 
the previous session we had to be content with upper and lower bounds. Although unmitigated fading 
reduces TC, it raises the possibility of opportunistic scheduling and transmit power control, which are 
discussed in §IV-DI and §IV-E[ 

A. General Fading 

With general fading values as in ([1]), the set of dominant interferers in ^ becomes 



Computation of the probability of a dominant interferer (P(n"(A) 7^ 0)) yields the following lower 
bound to OP ini: 



where the outer expectation is with respect to S. This expression is similar to the LB in (fT2l) . but 
the expectation in front of the exponential makes inverting this expression for A infeasible. Applying 



IV. Transmission Capacity in Fading Channels 




(30) 




(31) 
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Jensen's inequality to g (A) yields the following approximations: 

q{X) ^ 1 -exp|-A7rr2/3iE[/S]E[5~i]| (32) 

^ (33) 

7cr^(3^E[I^E[S-^ 

These quantities are approximations because Jensen's inequality yields inequality in the wrong direction. 
However, numerical results show that this approximation is reasonably accurate for small values of e 
[[T3l . It is possible to extend the upper bounds from §III-CI to fading [13J, but we focus exclusively on 
the above lower bound and approximation because they are more accurate. 

If we assume that the signal and interference coefficients follow the same distribution Fh, which is 
reasonable in most communication environments, the expressions in (|3T|)-(|33]) particularize to: 

g'(A) = I -Eh exp |-A7rr2/3-E[if^]i7-^ }] (34) 

g(A) ^ l-exp|-A7rr2/3^E[i/^]E[if-^]| (35) 

^ (l-e)log(l-e)-^ 
c(e) ^ 2 2 —■ (36) 

7ir^f3^E[H^E[H-^ 



Comparing the TC approximation in (1361) to the TC upper bound in (1131) we see that the effect of 
fading is captured by the term ^E[i7"]E[iJ^Q] j . By Jensen's inequality, this quantity is less than 
one (with equality only if H is deterministic) and thus fading has an overall negative effect relative to 
pure pathloss attenuation. Furthermore, note that the TC approximation in (|36l) is equal to the exact TC 
in (|39l) for Rayleigh fading derived in the next section. For the particular case of Rayleigh fading with 
a = 4, the approximate ratio (fT3l) over (l39l) equals | ^ 1.5708, while the exact ratio (([8]) over (l39l) ) is 
TT Q — ((i+e)/ 2)_^ which rapidly approaches | as e — )• 0. Thus, adding Rayleigh fading to a network with 



a = 4 reduces the TC by 57%. 
B. Rayleigh Fading 

The case of Rayleigh fading, where each Hij is exponentially distributed (unit mean), is appealing not 
only for its practical importance but also because it is one of the few cases for which the OP and TC can 
be computed in closed form. The following argument was made precise by Baccelli et al. [[24l. but can 
be traced to [|60l . [|6T1l . Define the aggregate interference seen at the origin as Z = X]ien(A) Hio\Xi\~°', 
and denote the Laplace transform of Z by Cz{s) = E [e^'*'^] . Then the success probability under 
Rayleigh fading is the Laplace transform of Z evaluated at s = /3r": 



POO 

P(SIR > /3) = F{Hoo > Pr'^Z) = / e-^""7z(2)d2 = E [e""^] 

Jo 



(37) 
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This transform can be computed explicitly, yielding an exact OP expression ((3.4) in n2M ) 



2 ^2 ItX I 27V 



g(A) = 1 -exp <j -A7rr^/3- — CSC ( — ) [> , (38) 



where esc denotes the cosecant. The corresponding exact TC expression is 



^ ^ 7rr2/^^^csc(^) ■ ^ ^ 



C. Nakagami Fading 



The Nakagami-m distribution has power given by 



, m m — 1 

m \ X / mx 



and is quite general in that Rayleigh fading corresponds to m = 1 and path loss only corresponds to 
m — 7- oo. Because the distribution is also of exponential form, OP and TC can be computed exactly in 
a manner similar to Rayleigh fading, resulting in a transmission capacity of [ITSl 



e 



1-1 



c(e) = — ■ 2 5 where (41) 



K, 



m—2 , k 



(42) 



(fc + 1)! J-1 

fc=0 ^ ' 1=0 

and B{a, h) = ^Y^J+h) ^el^ function. Although this expression is clearly more complex than (|39l ). 

it does describe nearly any fading environment. Interestingly, if m — )■ oo, i.e. for path loss only, PTI) 
converges to the upper bound of (fT3l) . 

D. Threshold scheduling 

Fading can potentially be exploited if only users experiencing good fading conditions transmit. This 
can be done through a simple threshold scheduling rule where each transmitter elects to transmit only 
if the signal fading coefficient i^oo is larger than a threshold t, as in lfT3l . Threshold scheduling is an 
example of opportunistic scheduling. The spatial intensity of attempted transmissions for threshold t is 
/i(t) = AP(ifoo > ^) = ^Fnit), i.e., the original intensity A thinned by the probability of being above 
the threshold. Because the threshold is on the received signal strength rather than the SIR, the decision 
depends only on local fading and does not affect the interference. Therefore, the outage probability with 
threshold t is: 

Hoo>t]. (44) 
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where the {Hij} are drawn iid according to Fh. The density of active transmissions is kept equal to 
u, independent of the value of t, by choosing A = Thus, the only change brought about is 

that the signal distribution follows distribution FH\H>t instead of Fh- As a result, the OP in (l44l) is 
decreasing in t and thus TC increases with The transmission capacity approximation is given by: 

. . (45) 

Comparing this with (l36l) . the (approximate) ratio of TC with threshold scheduling to that without it is 
— '^^'^2 — • Because bad signal fades are eliminated, the gains from threshold scheduling can be very 
substantial: for example, in Rayleigh fading a very reasonable threshold of t = 1 (i.e., dB) increases 
TC by a factor of 4.7, 3.3, and 2.25 for a = 2.5, 3, and 4, respectively. 



E. Power control 

While threshold scheduling attempts to completely avoid bad fades, an alternative strategy is to 
transmit regardless of the fading conditions and adjust transmit power to compensate for fading. In [fT4l 
a fractional power control policy in which each transmitter partially compensates for the signal fading 
coefficient is proposed. In particular, transmit power is chosen proportional to the fading coefficient 
raised to the exponent —7 where 7 G [0, 1]: 

ptx,fpc ^ P TT-J prx,fpc ^ P rrl-l -a 

Note that 7 = corresponds to constant power while 7 = 1 corresponds to full channel inversion. The 
resulting SIR is SIR = i^oV'^-"/ E.en{A) (^il'^io) |X,r". 

With channel inversion (7 = 1) there is no signal fading (S* = 1) and each interference coefficient is 
distributed as and thus based on (|3T1) we get the following OP lower bound: 

g''"(A) = 1 - exp |-A7rr2/3iE[/ft]E[if-^]| . (47) 

(There is no outer expectation because the signal fading coefficient is deterministic.) By Jensen's 
inequality, this quantity is larger than the OP lower bound for constant power given (|34l) . and thus 
the lower bounds indicate that inversion degrades performance. For Rayleigh fading this ordering is 
precise: the OP lower bound with channel inversion in (l47l) is equal to the actual OP with constant 
power given in (|39| ), and thus constant power is strictly superior to inversion in Rayleigh fading. 



'An outage is declared only if a transmitter actually attempts transmission and fails; not meeting the threshold is not considered an 
outage because it is essentially the same as not electing to transmit in pure Aloha. 
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Although inversion worsens performance, partial compensation for fading can be beneficial. If we 
consider general 7 and substitute the appropriate distributions for S and / in (l33l) . we get: 

,fpc. .... (l-e)log(l-e)-i 



Hi 



E 





E 






E 





(48) 

This approximation is maximized by minimizing E H~"'a E H^^^~"'^a over 7 G [0, 1], and an 
application of Holder's inequality yields 7* = 1/2. Although this only ensures that 7 = 1/2 is optimal 
for the TC approximation, results in [fT4ll confirm that 7 = 1/2 is also near-optimal for a wide range of 
reasonable parameter values Using 7 ^ | over-compensates for signal fading and leads to interference 
levels that are too high, while 7 ^ ^ leads to small interference levels but an under-compensation for 
signal fading. The benefit of FPC is substantial for small values of e and a. In Rayleigh fading, FPC 
increases TC by a factor of 2.1 and 1.2 for a = 2.5 and a = 4, respectively, for small e. 



V. Multiple antennas 

The amplitude and phase of fading channels vary quite rapidly over space, with an approximate 
decorrelation distance of half a wavelength (6 cm at 2.5 GHz). This allows multiple suitably- spaced 
antennas to be deployed at both the transmitter and receiver to generate NfNr Tx-Rx antenna pairs, 
where Nt and Nr are the number of transmit and receive antennas. Considerable work has been 
done on multi-antenna systems (MIMO) in the past decade, well summarized by [|62l|, ll63l . and such 
systems are now quite well understood and are central to all emerging high-data rate broadband wireless 
standards. However, much less is known regarding the use of antennas in ad hoc networks. In addition 
to providing diversity and spatial multiplexing benefits, multiple antennas also provide the ability to 
perform interference cancellation. Recent analysis of MIMO systems using the TC framework allows 
us to evaluate these different antenna techniques, and provides a very optimistic picture of the benefit 
of MIMO in ad hoc networks. 



A. Diversity 

Broadly defined, diversity techniques use TX and RX antennas to mitigate fading and increase the 
received SNR. With maximum-ratio combining/transmission (MRC & MRT), the transmitter and receiver 
apply weighting vectors at the antenna arrays based only on the Tx-Rx channel matrix. If the TX and 



^An important exception to this is for large values of e, i.e., dense networks, in which case the optimum tends towards constant power 
(7 = 0). 
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RX weight vectors are denoted by tg and tq, respectively, and Hj denotes the A^^ x Nt channel matrix 
from the i-th transmitter, then the SIR equation ([T]) becomes: 

SIR = . (49) 

Choosing the TX and RX weights as the right/left singular vectors of the largest singular value of Hq 
results in the signal coefficient being equal to the square of this singular value, and thus boosts signal 
power by a factor between max{A^t, A^^} and NtNr. With an appropriate application of (|33l) . this implies 
that the TC scales as ifBll : 

0{max{Nt, Nr}i) < c(e) < 0{{NtNr)i) as A^^, A^,, ^ oo. (50) 

The upper bound is tight for channels with high spatial correlation, while the lower bound is tight for 
i.i.d. Rayleigh fading. Note that A^^ = l,A^r > 1 and Nt > 1, Nr = I correspond to maximum-ratio 
combining (MRC) and maximum-ratio transmission (MRT), respectively. 

Orthogonal space-time block coding (OSTBC) is another diversity technique. OSTBC, which in- 
tuitively corresponds to repeating each information symbol from different antennas at different times, 
does not change the transmitted symbol rate but significantly increases received signal powerl^ However, 
interference power is also boosted and as a result OSTBCs increase the TC scaling only as c(e) = O^N" ) 
[[T5l . OSTBCs have very little affect on TC - the scaling gain is due to MRC at the receiver, independent 
of the code. 



B. Spatial Interference Cancellation 

If the receiver also has knowledge of the interferer channels, the A^r-dimensional RX weight vector 
can be used to cancel interference. In the single-transmit, multi-receive antenna setting with spatially 
uncorrelated Rayleigh fading, choosing the RX weight vector orthogonal to the vector channels of 
the strongest A^^. — 1 interferers (i.e., tq -L Hi, . . . , HAr^_i) results in ©(A^^^ °) TC scaling [fT6l . 
An even larger TC increase is obtained if the RX vector is designed to cancel interference and reap 
diversity. In particular, using about half the RX degrees of freedom for cancellation and the remainder 
for diversity (i.e., choosing vq as the projection of vector Hq on the nuUspace of Hi, . . . , HAr,,/2 ) 
leads to 0{Nr) TC scaling IfTOl B In fact, the SIR is maximized, and thus the benefits of interference 

^For some combinations of Nt and Nr OSTBCs either lose orthogonality, or reduce the data rate slightly. The results here make the 
optimistic assumption of rate 1 orthogonal STBCs for general Nt,Nr. 

■*Both of these scaling results are obtained using the OP upper bounding techniques described in i]III-CI 
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cancellation and diversity are optimally balanced, if the RX vector is chosen according to the MMSE- 
criterion: = ( 2Jjen(A) °HjH] j Hq. The MMSE filter is generally quite difficult to deal with 
analytically, although large-system results are derived using random matrix theory in ||64l . 

In Fig. |4] the TC of diversity (beamforming and OSTBC) and intererence cancellation are plotted 
versus the number of antennas (A^) for a = 4 and /3 = 1. All of the techniques except OSTBC provide 
significant gains, but the combination of interference cancellation and diversity clearly provides the 
largest TC, as predicted by the TC scaling results. 

C. Spatial Multiplexing 

The most aggressive use of the antennas is to use them to form up to L < mm{Nt, Nr} parallel spatial 
channels. If the transmitter has knowledge of the channel matrix Hq, this corresponds to beamforming 
along the eigenmodes of the channel. The achieved SINR for each spatial channel depends on the 
eigenvalues of the channel matrix as well as the interference power, so some channels are much better 
than others. When subject to an SINR target and an outage constraint, it is preferable to transmit only a 
small number of streams (L <^ N) unless the network is very sparse. This is illustrated in Fig. [5] where 
the optimized number of spatial streams (as determined in GUI ) is plotted versus the interferer density 
and this quantity is seen to decrease from to 1 with the density. Ideally, the number of spatial channels 
can be adapted dynamically based on the channel and interference strengths to maximize the quantity 
Lc{e,L), which is the area spectral efficiency (ASE) shown in Fig. [51 and has a unique maximum 
EOl . Here c(e,L) is the TC with target OP e when L antennas are employed. If each TX wishes to 
communicate with multiple receviers, multi-user MIMO techniques can be used to send separate data 
streams to each receiver. In the situation where each transmitter and receiver has N antennas, the TC has 
been shown to increase super-linearly with when dirty paper coding, the optimal multi-user MIMO 
technique, is used [TSl- 

If the transmitter does not know channel matrix Hq, spatial multiplexing is generally performed by 
transmitting independent data streams from each transmit antenna. The OP and TC for low-complexity 
(and sub-optimal) MRC and zero-forcing receivers are known [8], but many important questions remain 
unanswered on this topic, e.g., performance with optimal MIMO receivers. 

VI. Current Limitations and Future Directions 

Although the results presented in this paper have illustrated the value of the transmission capacity 
framework, they have also failed to capture two important aspects of ad hoc networks. The first is 
that they are for a snapshot, or single-hop, of the network. This may be acceptable for unlicensed 
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spectrum analysis or other decentralized networks, but ad hoc networks must route traffic from source 
to destination, often over multiple hops through intermediate nodes. A network with higher single-hop 
TC should be able to achieve higher end-to-end capacity than a network with smaller TC because more 
simultaneous transmissions are possible. However, important issues such as desired hop length, number 
of hops, multi-hop routes, and end-to-end delay are not presently addressed. In addition, noise should 
not be neglected since a principle function of multihop is to increase the SNR for each hop. Some work 
that attempts to use the results of this paper (or similar results) to address multihop includes [|24ll . where 
a metric called expected forward progress is introduced and used to find the optimum split between 
transmitters and receivers (potential relays) in terms of the Aloha contention probability. Recently, ll65l 
has developed a multihop model and found an end-to-end delay-optimizing strategy in a Poisson field 
of interference (without noise), while fi66l1 finds the end-to-end transmission capacity in closed-form 
(i.e., transport capacity) with noise under a few restrictive assumptions like equi-distant relays and 
independent retransmissions. Clearly, this is a line of work that should be pursued and improved upon 
in the coming years. 

The second lacking aspect of the current results is that they rely on a homogeneous Poisson distribution 
of nodes for tractability, which accurately models only uncoordinated transmissions {e.g.. Aloha). 
A well known alternative is to schedule simultaneous transmissions with the objective of control- 
ling interference levels. Local scheduling mechanisms generally space out simultaneous transmissions, 
thereby significantly changing the interference distribution, while idealized centralized scheduling can 
eliminate outages altogether and determine the optimal set of transmitters in each slot {e.g., max- 
weight scheduling within the backpressure paradigm [|67l). Preliminary work in this direction includes 
computing the outage probability and transmission capacity under non-Poisson point processes [|4T| . [|68l . 
[|69l . Although scheduling mechanisms provide obvious gains, these come at the cost of overhead {e.g., 
control messages). Thus, a general open question is understanding the tradeoff between the benefits 
and overhead costs of different scheduling/routing mechanisms (Aloha is a particular point on this 
tradeoff curve), and determining the appropriate techniques for different network settings. Furthermore, 
a fundamental property that applies even to scheduled systems is that transmissions occupy space 
whenever interference is treated as noise; the transmission capacity provides a clean characterization of 
this space, and thus many of the insights apply, in principle, to scheduled systems as well. 

As is true of any complicated research topic, discussion of a particular model or framework exposes 
tension between analytical tractability and accuracy/generality. The transmission capacity framework 
clearly leans towards simplicity and tractability, but nonetheless provides valuable design insight and a 
launching point for more refined, less tractable network analyses. 
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TABLE I 

Notation used in paper. 



a = b 


a is defined to equal b 


A 


spatial intensity of attempted transmissions (per m^) 


n = {X.} 


Poisson point process (PPP) of intensity A of transmitter locations 


a 


pathloss exponent (a > 2) 




SIR/SINR requirement for successful reception 


r 


distance separating each Tx-Rx pair 


iW 


outage probability (OP) 


e 


constraint on OP 


c{e) 


transmission capacity (TC) 


P 


transmission power 


Hij 


fading coefficient from transmitter i to receiver j 


M 


number of frequency channels, or spreading factor 


Nr,Nt,N 


number of receive, transmit, or total antennas 




Fig. 1. The transmitter locations (black circles) at a typical time form a Poisson process, II; each transmitter has an assigned receiver 
(gray circles) located at distance r. The reference communications link has a reference receiver at the origin (green) and a reference 
transmitter at distance r (red). Each black transmitter generates interference seen at the reference receiver, indicated by the dashed lines. 
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q(A) 




Fig. 2. Top: OP g(A) versus the spatial intensity of attempted transmissions, A, for tlie basic model with a = 4, P — 3, and r — 10 
meters. The three lines are lower bound, exact OP, and the (Chemoff) upper bound. Bottom: the TC c(e) versus the outage requirement 
e obtained by inverting the outage expression and bounds. 
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Fig. 3. Top: Optimized SINR threshold P versus path loss exponent a. Bottom: Optimized outage probability constraint e versus path 
loss exponent a. 
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Fig. 4. The transmission capacity of various spatial diversity techniques vs. the number of antennas per node. Interference cancellation 
& diversity refers to cancelling the nearest N (2 interferers and using the remaining degrees of freedom for diversity. Here, e = .1, a = 4, 
^ = 1. 




Fig. 5. Optimal number of MIMO modes L and Area Spectral Efficiency (ASE) vs. transmitter intensity per vc? . The L curves are 
monotonically decreasing, ASE curves are bell-shaped and have a unique maximum. Here, e = .1, r = Im, a = 4. 



